from MATRIX_INV_DET import matrix_inv,matrix_det
import numpy as np
import scipy as sc
# Test matrix
A = np.array([[1,2,3],[0,1,4],[5,6,0]], dtype=np.float64)

# Calculate inverse and determinant
A_inv = matrix_inv(A)
A_det = matrix_det(A)

print("Inverse of A:")
print(A_inv)
print("Determinant of A:")
print(A_det)

print(f"numpy.linalg.inv(A):\n{sc.linalg.inv(A)}")
print(f"numpy.linalg.det(A):{sc.linalg.det(A)}")